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Abstract. We report a theoretical study of the many-body effects of electron- 
electron interaction on the ground-state and spectral properties of bilayer 
graphene. Using a projector-based renormalization method we show that if a 
finite voltage difference is applied between the graphene layers electron-hole pairs 
can be formed and — at very low temperatures — an excitonic instability might 
emerge in a bilayer graphene structure. The single-particle spectral function 
near the Fermi surface exhibits a prominent quasi-particle peak, different from 
neutral (undoped) bilayer graphene. Away from the Fermi surface, we find that 
states strongly interact with plasmons, thereby giving rise to a broad plasmaron 
peak that explains prominent features of recent angle-resolved photoemission 
spectroscopy data. 



PACS numbers: 71.10.-w, 71.38.-k, 75.47.Gk, 71.70.Ej 

1. Motivation 

Graphene based structures are most likely the building blocks of future nanoelectronic 
devices. The reasons for this are manifold, but one may highlight that many of the 
exceptional properties of this new class of low-dimensional materials, which arise from 
the special form of energy spectrum near the so-called Dirac nodal points and the 
related nontrivial topological structure of the wave function, can be easily modified 
by the application of external electric and magnetic fields, as well as by confining 
the sample geometry, by chemical doping, or by edge functionalization and substrate 
manipulation; for recent reviews see p~H3]. 

Single layer graphene (SLG) [IJ[5], a truly two-dimensional (2D) crystal with 
remarkable mechanical properties, can be considered as a gapless semiconductor with 
zero density of states at the Fermi level and a linear energy dispersion close to the 
(inequivalent) corners of the Brillouin zone (K, K' Dirac points). Thus the low-energy 
electrons are massless, chiral Dirac fermions. As a consequence any backscattering is 
suppressed and the charge carriers are almost insensitive to disorder and electron- 
electron interactions. Bilayer graphene (BLG) [HE], consisting of two graphene 
monolayers, is also a zero-gap semiconductor, if unbiased, but with a parabolic 
band dispersion. That is, in this material the low-energy electrons acquire a finite 
quasiparticle mass. 
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In both SLG and BLG, the charge carrier density and polarity can be controlled 
by simple application of a gate voltage |5]. This electric field effect is fundamental 
for potential technological applications. For BLG even the band structure might be 
manipulated by an electric field, so that a gap between the valence and conduction 
bands varies between zero and mid-infrared energies [9]. 

BLG, besides being at present the only known semiconductor with a tunable band 
gap, is intriguing also from a many-particle interaction physics perspective. While in 
SLG the interaction parameter is the fine-structure constant and Coulomb interaction 
effects are typically small, in BLG the interaction strength is alterable by changing 
the carrier density [10]. The ability to modify the carrier polaribility of an individual 
layer even implies that the interlayer Coulomb interaction can be converted from 
repulsive to attractive [11]. Since a gate bias across the two-layered graphene leads 
to a charge imbalance in the layers, an attractive Coulomb interaction between the 
excess electrons and holes on opposite layers raises the possibility of formation of 
electron- hole bound states (excitons) [12] . Excitonic effects have been reported in the 
optical response of doped SLG [13] and BLG [13] . Interestingly, BLG seems to be 
also a promising candidate for electron- hole pair condensation [TTj[T5UT7] . In the weak 
coupling regime, exciton condensation is triggered by a Cooper-type instability of 
systems with occupied conduction-band states and empty valence-band states inside 
identical Fermi surfaces p~5j[T8]. The particle- hole symmetry of the Dirac equation 
obviously ensures a perfect nesting between the electron Fermi surface and its hole 
counterpart in the n- and p-type layer of a biased BLG system [15 . A condensate of 
spatially separated electron-hole pairs then might occur when the tunneling between 
the layers is negligible, but interlayer Coulomb interaction is not [TT| [T6]. Such bilayer 
exciton condensates typify as counterflow superfluids [T9] . 

In characterizing many-body aspects of BLG, and particularly the electron- 
electron, electron-hole and electron-plasmon interaction effects, the single-particle 
spectral function A(k,u) gives detailed information. Angle-resolved photoemission 
spectroscopy (ARPES) is a powerful probe of A(k,u) in (quasi-) 2D single-, bi- 
and few-layer graphene because it achieves frequency and momentum resolution and 
thereby directly addresses the quasiparticle energy band properties [jJ|20j[2T]. So far 
the spectral function of SLG and BLG has been calculated within the simple GqW 
and random-phase approximation [10.22-25 ], with a focus on dynamic screening, Kohn 
anomaly, Friedel oscillations and plasmons. 

In this work, we employ the projector-based renormalization method (PRM) to 
determine the ground-state and spectral properties of an effective graphene bilayer 
model describing the interaction and possibly pairing of electrons from the top bilayer 
with holes from the bottom layer [cf. figure [T] (a)]. Our main interests here are 
the effects of the Coulomb interaction regarding the formation of excitonic bound 
states and, for the doped system, the coupling between the charge carriers and 
plasmons leading to composite "plasmaron" quasiparticles. We show that graphene- 
based bilayers might allow for electron- hole pair condensation and exhibit, in their 
ARPES spectra, pronounced plasmaron signatures away from the Dirac point. 

2. Theoretical Approach 

We consider two graphene sheets separated by a dielectric thickness and assume a 
hexagonal stacking in which each sublattice in one layer is on top of the corresponding 
sublattice in the other layer [TT] , Let now V g = eE e ^ t d be a gate- induced potential 
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Figure 1. (a) Schematic of a graphene bilayer system with electron and hole 
carriers induced by external gates forming excitonic bound states due to their 
Coulomb attraction, (b) The positions of the chemical potential in the two 
graphene layers is adjusted by the voltage V g . (c) Brillouin zone of the graphene 
honeycomb lattice (black lines) and equivalent Brillouin zone (solid green lines) 
spanned by the inverse lattice vectors Gi and G2. Dirac cones are located at the 
K and K' points. 



difference between the two layers. Then the external electric field £? ext outside the 
bilayer system can be used to tune the chemical potential fi = V g /2 [15] : For a neutral 
graphene bilayer system the Fermi level lies symmetrically in the conduction band 
of the upper (n-type) layer and in the valence band of the lower (p-type) layer [see 
figure [D (b)]. Within an effective two-band description, the filled valence and empty 
conduction bands of the n- respectively p-type layers can be neglected [dashed lines in 
panel (b)], and the Hamiltonian of the graphene bilayer, in the absence of the Coulomb 
interaction, is 

^0 = 9s e i a k a k + 9s £ k 6 k 6 k • (!) 

k k 

Here and are the creation operators of electron and hole quasiparticles with 
in-plane momentum k and band dispersions 

1/2 



±7o 



1 + 4 cos — + 4 cos — cos — k x 
2 2 2 



(2) 



which for momenta close to the Dirac points become the linear energy-momentum 
relations sketched by the blue and red straight lines in figure [T] (b), respectively. In 
equations |I])-(j2j), g s = 2 allows for the spin degeneracy and 70 denotes the nearest- 
neighbor transfer amplitude. 

In order to describe exciton formation, we take into account the interlayer coupling 
between electrons and holes, 



k 2 q ^ki+q^ki^ks-q^ks 



kik 2 ,q/0 



where 



kik 2 q 



e" rf|q| 01 4>2 ... 
v — j — j — cos — cos — with k = g s 



2ire 2 



(3) 



(4) 



The Coulomb matrix element /7k 1 k 2 q contains besides the dielectric constant e, 
characterizing the embedding medium, a graphene-specific factor oc cos 4r cos 4r , with 
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scattering angles fa = 0^ — 6^+q where 0^ = a,taxi(k y / k x ) [11 . At this point at least 
two remarks are in order. First, we have omitted the intralayer Coulomb repulsion 
between electrons or holes that manifests itself simply by a screening of the interlayer 
Coulomb interaction. For bilayer graphene it has been shown that the dynamically 
screened interlayer electron- hole interaction is attractive [IB] . Second the effective 
Coulomb attraction is long-ranged and depends in an involved way on the transmitted 
momentum. This is a major difference to the momentum-independent short-ranged 
Coulomb attraction in the extended Falicov-Kimball model, recently studied in the 
context of exciton formation and condensation within the PRM approach [26] [27]. 
Excluding the q = component (comprising the jellium background), the Hartree 
contributions to the one-particle energies vanish. 

To proceed, we rewrite the total Hamiltonian, H = Ho + ^e-p? m a normal- 
ordered form 

n = s H a ^ k + + ( A ^k% + h.c.) 

k 

+ Jf ^kaq : 4 1 +q ft k 1 ^ 2 -q 6 k 2 (5) 

kik 2 ,q/0 

and look for a non- vanishing excitonic expectation value 

k e~ d ^ (1 + cos0) . | 7 . 
A * = ~n E ~\^~ 2 <4 +q W > (6) 

indicating a spontaneous symmetry breaking due to the pairing of electrons and holes. 
T-L is given in units of g s , and <p = #k+ q — #k- 

Now the single-particle spectral function, e.g. for the conduction band electrons, 

A a (k,uj) = --lmG a (k,uj) (7) 

7T 

with G a (k, uo) = ((^k' a k))^ being the Fourier transform of the retarded Green 
function (uj — >> cj + i0 + ), can be computed by the PRM in close analogy to the 
calculation performed for the extended Falicov-Kimball model [26] . We have 

G a (k,co) = ((a k ;al)) fL , (8) 

where H is a renormalized Hamiltonian with all transition energies due to the 
interaction part of the Hamiltonian (|5]) are successively integrated out. After 
this renormalization procedure becomes complete, H embodies a noninteracting 
Hamiltonian, but with a renormalized band structure, and therefore can be easily 
diagonalized by a Bogoliubov transformation to yield 

H = J2(KaW + Elblb k ) , (9) 



where 

K / b = _ /+ sgn( £ ^- g+) Wk+ e£+K with Wk= [ (e -+_ e -- ) 2 + 4 | Ak |2j 1/2 . (10) 



The new fermionic operators satisfy a k = ^ k a k + ^k^k anc ^ ^k = v k a k — ^k^k> with 
^k + v k = 1j v k = — s S n (^k — ^k)(^k — ^k)/^k]- Replacing the renormalized 
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quasiparticle operators = x^a^ + S pq ^kpq a k+q : ^p-cfip : m equation (|5J), the 
spectral function is given to leading order in A by 

A°(k, u) = \x k \ 2 [vlS(w - E£) + i£5(u - El)} + ^ £ l^ Pq | 2 x (H) 

pq 

xS(u- (££ +q - ^ + £*_ q )) [<n£ +q >«n» _ q > - <n*» + <n£)(l - <n*_ q ))] . 
Note that the expectation values 

(n£)=4/(^) +u 2 /( ^) a nd (4) = t£/(23£) + (12) 

have to be calculated with the fully renormalized Hamiltonian H as well, where f(E^ b ) 
are Fermi functions. The spectral function for the valence-band holes, A 6 (k, u), is 
obtained in the exact same manner. 

3. Results and Discussion 

The PRM approach outlined in the previous section must be numerically evaluated 
of course. In doing so, we work in the reciprocal space bounded by the green lines 
in figured] (c), on a discrete set of N = 64 x 64 points. Convergence of the whole 
renormalization scheme is assumed to be achieved if all quantities are determined with 
a relative error less than 10 -5 . Thereafter the spectral functions are calculated using 
a Gaussian broadening in energy space of width 0.06. In what follows, we choose 
7o = 1 as unit of the energy, measure all energies relative to the Dirac-point chemical 
potential of the balanced bilayer, and fix the inter layer Coulomb interaction strength 
to be k = 7. 
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Figure 2. Renormalized band dispersions, E^ ,b , along the high-symmetry 
directions of the Brillouin zone of graphene, indicated in figure^ by the dashed 
green lines. Results are given for an external field E eyi t = 1 at zero temperature 
(solid lines) and at T = 0.1 > Tsf (dashed lines). 

Figure [2] shows how the band structure of a suchlike interacting BLG system 
changes, if the electric field £^ e xt = 1 is applied. Apparently a pronounced excitonic 
gap appears at zero (low) temperatures Q The gap originates from electron- hole pair 



X Note that the natural stacking in graphite is the Bernal stacking which merely leads to a 
insignificantly smaller excitonic gap however, see figure 2 in 
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(exciton) condensation, driven by the attractive Coulomb interaction over the full 
bandwidth. In the result of this excitonic instability a macroscopically coherent 
quantum state emerges which, since electrons and holes are spatially separated in 
a graphene bilayer, differs in nature from the excitonic insulator phase in bulk 2D and 
3D systems [27H35] . As can be seen from figure [2j the renormalized quasiparticle bands 
E^ b possess a nearly perfect particle-hole symmetry. The gap disappears at higher 
temperatures when inter layer quantum coherence is destroyed by thermal fluctuations. 
For a strictly 2D system, of course, the critical temperature for exciton condensation 
would be zero, but the superfluid properties should survive for temperatures smaller 
than the Kosterlitz-Thouless transition temperature [Tl][T5] . Therefore, for the BLG 
system, our finite-temperature results will be valid for T < Tkt at least. 

The zero- (low-) temperature exciton condensate is characterized by a non- 
vanishing expectation value 5^ = (a^.& k ). Since the interlayer Coulomb interaction 
between electrons and holes is attractive for all momenta and we exclusively consider a 
uniform exciton condensate in BLG, S = (l/N) Xlk( a k^k) can ^ e taken as an adequate 
order parameter for the electron- hole pairing. This quantity is plotted in figure [3] (a) 
against the magnitude of the electric field at zero temperature (red squares). It shows 




Figure 3. (a) Order parameter S = (l/N) X]k( a k^k) as a function of temperature 
T for an external electric field E'ext = 1 (black circles) respectively as a function 
of Eext at T = (red squares). The other panels show the magnitude of the 
expectation value, 5^ = (o^b^) with k = /cj.Gi + /C2G2, in reciprocal space by an 
intensity plot for E ex t = 0.6 , T = [panel (b)], E ex t = 1, T = [panel (c)], and 
Eext = 1, T = 0.1 [panel (d); note the different intensity-color coding]. 
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that above a certain threshold the excitonic condensate becomes more robust as more 
charge carriers (electrons and holes) are induced into the layers by increasing the bias 
between them. Obviously, a larger valence and conduction band overlap amplifies the 
Cooper- like excitonic instability in the perfectly nested BLG system. Fixing E ext = 1, 
the order parameter is suppressed as the temperature increases (black circles), and 
vanishes when the electron-hole coherence is completely lost in the high-temperature 
phase. Panels (b)-(d) yield detailed information on the variation of the the order 
parameter function 5^ in the BLG Brillouin zone (here the diagonal belongs to the TM 
direction). Raising the external electric field at zero temperature, exciton formation 
and condensation can set in if the conduction band of the upper layer enters the 
valence band of the lower layer. This takes place first at the Dirac points. If the 
electric field is enhanced further, electron and hole bands strongly hybridize along a 
crater-rim shaped contour around K, K', thereby establishing electron- hole coherence 
in this region [see panels (b) and (c)]. Above a critical temperature T > T$ F this 
signature is washed out and 5^ becomes basically zero anywhere. 

To gain further insight about the Coulomb interaction effects in BLG we now 
analyze the PRM results obtained for the total wave-vector and frequency resolved 
single-particle spectral function A(k,oj) = A a (k, uj) + A 6 (k, uj). Depending on 
the energy a;, A(k,uj) measures the probability of finding an electron (hole) with 
momentum k in the system and is directly accessible by photoemission (inverse 
photoemission) experiments. Thereby it gives valuable information about the 
quasiparticle properties, particulary with regard to the quasiparticle dynamics. 

A key question in this respect is whether the non- Fermi-liquid behavior observed 
for neutral undoped BLG at the Dirac point will survive in the presence of any doping. 
While Coulomb interactions self-evidently play no role for a pure zero-density Fermi 
system, i.e. will not affect this non-Fermi-liquid fixed point, they become important 
at any finite carrier density (induced, e.g., by applying a gate voltage), or if spatially 
density fluctuations are present in real BLG samples. Using the random-phase and 
parabolic band approximations, it has been shown quite recently, that BLG constitutes 
a Fermi liquid for doped cases, similar to the more standard 2D semiconductor-based 
electron gas [TO] , If so, the Coulomb interaction should cause noticeable correlation 
effects. 

Figure [4] gives the intensity of the single-particle spectral function, subject to uj 
along the high- symmetry directions of the Brillouin zone, for several characteristic 
situations. Panel (a) displays the behavior of A(k, uj) when a exciton condensate is 
realized at T = 0. We find the location of the quasiparticle bands in accordance with 
the results of figure [21 These signatures have large spectral weight. Accordingly the 
effective-mass renormalization in BLG is expected to be rather weak, in contrast to 
the large mass renormalization in a normal 2D electron gas [10]. The single-particle 
excitation gap, located in the vicinity of the K point, disappears for T = 0.1 [panel (b)]. 
That is, for T > T$f we have gapless excitations only. Panels (c) and (d) show how 
the overlap of the predominantly electron (uj > 0) and hole (uj < 0) quasiparticle 
bands shrinks as the electric field is reduced at fixed temperature (T = 0.1). 

Besides the quasiparticle band signal discussed so far, we observe — especially in 
panels (a)-(c) — another signature that has less spectral weight, but is nevertheless 
well pronounced. It can be attributed to a plasmaron [36] . In an electron gas, the 
long-range Coulomb interaction manifests itself through plasmonic collective charge 
density oscillations. These plasmons can in turn interact with the charge carriers, with 
the result that a new composite particle called plasmaron is formed [36]. Plasmaron 
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quasiparticles have been predicted to occur in SLG [22,37,38] and BLG [10], and 
recently were observed in the former system by ARPES [21 j. Our, regarding the 
treatment of interaction and lattice structure effects more involved PRM calculation 
corroborates this prediction for graphene bilayers. As can be seen from figure [H 
satellite bands, arising from the dressing of the injected particle (electron or hole) 
by plasmon modes, appear away from the Fermi momentum. The plasmaron shows 
up at greater (smaller) frequency than the quasiparticle because of the extra energy 
cost for the binding of the electron (hole) with the plasmon. The higher intensity 
of the plasmaron in the lower branch of the spectrum can be attributed to the fact 
that the plasmon oscillations are more pronounced in the p-type (few- hole) layer. The 
plasmaron signature weakens, of course, at very high doping level. Increasing the 
temperature, the plasmaron signal is slightly enhanced initially (since the electron- 
hole pairing is suppressed), but finally is smeared out when the temperature achieves 
the order of magnitude of the particle-plasmon binding energy. For small gate voltages 
[panel (d)], the plasmaron dispersion resembles that of SLG 137. 38^. simply because 
the band overlap is small or even absent and the quasiparticles in the n- and p-type 
branches become uncorrelated. 
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Figure 4. Intensity plots of the total wave- vector resolved single-particle spectral 
function A(k, oo) for T = and E e *t = 1 [panel (a)], respectively T = 0.1 and 
£ext = 1 (b), 0.6 (c), and 0.2 (d). 
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Figure [5] finally presents the frequency dependence of the partial electron and 
hole spectral functions at selected k points. It shows that the quasiparticles far away 
from the Fermi point are predominantly hole-like (red curves) or electron-like (blue 
curves) at both zero and finite temperature; cf. the (a) panels. This holds as well 
as for the corresponding plasmaron signatures. For k ~ kp (Fermi wave vector), the 
single-particle excitations are gapful (gapless) for T = (T > Tsf)- Most notably, 
we observe a strong mixing of electron and hole bands, which also becomes apparent 
for the plasmarons. Moving to the K-point, this mixing is lost [cf. (c)-panels], and 
the plasmaron bump appears well below (above) the corresponding electron (hole) 
quasiparticle peak. 
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Figure 5. Single-particle spectral functions of conduction-band electrons 
A a (k, uo) (blue lines) and valence-band holes A b (k,w) (red lines) plotted against 
uj for characteristic momenta k = Y [left-hand (a)-panels], k ~ kp [middle (b)- 
panels], and k = K [right-hand (c)-panels] at T = 0.0, E e ^t = 1.0 (upper row 1) 
and T = 0.1, £ ex t = 1.0 (lower row 2). 



4. Concluding remarks 

Analyzing by means of a projector-based renormalization method Coulomb interaction 
effects within a generic graphene bilayer model, we suppose that electron-hole pair 
formation and belike condensation might appear in a BLG system at least at zero 
temperature, if charge imbalance between the layers is induced by an external electric 
field. Thereby the condensation of the spatially separated, but bounded electrons and 
holes turns out to be a robust phenomenon that will not require a precise alignment 
of the graphene sheets embedded in the dielectric [TTJ[T5] • The excitonic instability is 
of Cooper (BCS) type, and manifests itself by an arising interlayer coherence between 
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the electron and hole quasiparticles residing inside nearly perfectly nested Fermi 
surfaces. The condensed superfluid low-temperature state in BLG differs from the 
excitonic insulator state observed in (intermediate valent) rare-earth [30, or transition- 
metal systems [33]. The exciting BCS-BEC crossover upon an increase of the 
coupling strength [3T j l34 l l35] , discussed for the latter materials as well as for standard 
semiconductor-based coupled quantum well structures, seems to be achievable in BLG 
in the presence of a perpendicular magnetic field only, when localized magnetoexcitons 
will form (see [16] and references therein). The pronounced quasiparticle peak in the 
single-particle spectra indicates that the electrons (holes) in doped BLG behave more 
like in a usual Fermi liquid. The ARPES spectra, in addition, feature a plasmaron 
peak, arising from the coupling of the charge carriers to collective density oscillations. 
This effect might be of importance for the functionality of new plasmonic devices 
merging photonics and electronics [39) 140]. 
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